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The energetics of thin elongated ferromagnetic nano-islands is considered for some different shapes, 
aspect ratios, and applied magnetic field directions. These nano-island particles are important for 
artificial spin-ice materials. For low temperature, the magnetic internal energy of an individual 
particle is evaluated numerically as a function of the direction of a particle's net magnetization. 
This leads to estimations of effective anisotropy constants for (1) the easy axis along the particle's 
long direction, and (2) the hard axis along the particle's thin direction. A spin relaxation algorithm 
together with fast Fourier transform for the demagnetization field is used to solve the micromagnetics 
problem for a thin system. The magnetic hysteresis is also found. The results indicate some 
possibilities for controlling the equilibrium and dynamics in spin-ice materials by using different 
island geometries. 
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I. INTRODUCTION: ELONGATED THIN 
FERROMAGNETIC NANO-ISLANDS 

Disordered and frustrated magnetic states such as 
those present in artificial spin ices^ continue to draw 
interest, due to their competing ground states, magnetic 
monopole excitations^, string excitations^— and the diffi- 
culty to achieve thermal equilibrium. Those systems are 
composed from elongated magnetic islands or particles of 
some length L x (several hundred nanometers) and width 
L y grown or etched lithographically to a small height L z 
on a substrate, whose geometric demagnetization effects 
(effectively, internal dipolar interactions) lead to a strong 
magnetic anisotropy. The typical islands have L z much 
less than L x or L y . Obviously any very thin magnet 
acquires an effective easy-plane anisotropy^, and if the 
particle is narrow as well, the long direction becomes an 
easy axis. The demagnetization field within an individ- 
ual particle is responsible for this, making the plane of 
the island (xy-plane) an easy plane, and the x-axis an 
easy axis. Then net magnetic moment fx acts somewhat 
like an Ising variable with a defined easy axis x. These 
islands are arranged into ordered arrays to produce, for 
example, square lattice or kagome lattice artificial spin- 
ices. The analysis of spin ice models assumes that such 
particles have only the two states with fx either aligned or 
anti-aligned to the particle's easy axis. The dipolar inter- 
action between different particles on one of the spin-ice 
lattices leads to the ice-rules, such as the "two in / two 
out" rule for the square lattice and pyrochlore spin ices^. 
Such ice rules are only energetic preferences, however, 
and only indicate the preferred states of the magnetic 
moments. They are not absolute rigid statements about 
the allowed states. Thus, the intention here is to inves- 
tigate the energetics of the fluctuations away from this 



Ising aligned state, in the individual elliptical islands that 
are used to compose a spin-ice system. 

At some level, there must be transitions between these 
Ising-like states. An individual particle may contain 
thousands of atomic spins, leading to a substantial en- 
ergy barrier that must be surpassed to flip the Ising state 
of a particle. Hence, the dynamics is greatly constrained 
by such energy barriers. It is our interest here to discuss 
how this barrier depends on the particular geometry of 
the islands, and make some evaluations of the depen- 
dence of the effective potential on the island shape and 
height. The types of shapes we consider are ellipses. Thin 
single domain ellipses were studied by Wei et al£, who 
found that the reversal process involves close to a uniform 
Stoner-Wohlfarth rotation, but with reduced energy bar- 
riers due to some non-uniformity of the magnetization. 
However, we find here that for high- aspect ratio ellipses, 
this non-uniformity is minimal and a uniform rotation 
model could be very useful. 

Although the theory for spin ice has been developed for 
Ising-like magnetic moments, their dynamics requires a 
different model. In reality, the underlying magnetic mo- 
ment must be evolving from much more complex dynam- 
ics. The reversal of an individual island, in the dipolar 
fields of its surrounding islands, must be a complex pro- 
cess, and could involve the motion of domain walls and 
vortices within the individual particles, or an impeded ro- 
tation of the local magnetization mostly in unison. But in 
the assumption of strong ferromagnetic exchange inside 
a particular particle, and a uniform externally applied 
field, one can investigate the reversal process using differ- 
ent approaches to the micromagnetic&i^, and see whether 
vortices or domain walls play any significant role. Espe- 
cially, one can investigate whether there are intermediate 
metastable vortex or domain-wall states as steps of the 
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reversal. To a great extent for the thin elliptical particles 
considered here, the reversal proceeds mostly as a nearly 
uniform but impeded rotation of the magnetization of 
the particle^, although the switching fields are reduced 
compared to a perfectly uniform rotation. Hence, the 
idea of an Ising spin for a particle can be replaced by a 
three-dimensional magnetic moment ft, moving in some 
anisotropy potential, but free to point in any direction, 
if enough energy becomes available to it. 

Obviously, by changing the aspect ratios g\ = L x /L y 
and g 3 = L x / L z of the particle, its effective anisotropy 
changes. The deviation of the ratio L x / L v from 1 deter- 
mines the strength of an easy-axis anisotropy constant, 
call it K\ , for the net magnetic moment to rotate within 
the xy plane. The other aspect ratio of length to thick- 
ness, L x /L Zl determines the difficulty for the magnetic 
moment fl to tilt out of the xy plane. Thus it deter- 
mines the strength of a hard axis anisotropy constant, 
call it K3. The goal here is to make some accurate esti- 
mates for these constants, and in the process, to justify a 
more generalized description of the magnetic dynamics, 
not based on an Ising variable, but rather, on an effective 
three-dimensional magnetic moment, that is allowed to 
make deviations from the Ising axis. For a particle whose 
hard axis is along z and easy axis is along i, an effective 
potential that approximately represents their energies is 
shown to be 

E = E + K 1 [l-{fi-x) 2 ]+K 3 (fi-z) 2 (1) 

where fi is the unit vector pointing in the direction of the 
particle's net magnetic moment. Eq is the energy when 
the magnetic moment fl is along the easy axis. This type 
of potential is continuous, in contrast to the two-state 
Ising particle, having a well-defined energy barrier, and 
having a more realistic dynamics. Further, it will give the 
possibility of controlling the thermodynamics of spin-ices 
via changes or variations in the nano-island structure, 
that can modify the energy barrier. 

The calculational approach is a modification of usual 
micromagnetic o 1 1 ' 1 2 , as follows. A particle is partitioned 
into cells of size a x a x L z , under the assumption of 
the local magnetization M(r) being independent of the 
^-coordinate (along the thin dimension). Thus, there is 
only a single layer of cells in the a;y-plane, with the de- 
sired shape, say, an ellipse of major diameter L x and 
minor diameter L y < L x . The saturated magnetization 
in each cell interacts with the neighboring cells by ferro- 
magnetic exchange, an externally applied magnetic field, 
and interacts with all cells via the demagnetization field. 
The demagnetization field is calculated using an effective 
Green's function that applies for thin systems^, see be- 
low, with the calculation accelerated by using a 2D fast 
Fourier transform (FFT). To evolve towards the nearest 
(possibly meta-) stable magnetic state, we do not use 
integration of the Landau-Gilbert spin dynamics equa- 
tions with damping. Instead, a faster procedure is to 
use a local spin-alignment algorithm, that involves no 
damping parameter. In one step of this algorithm, each 



cell's magnetic moment is pointed towards the local total 
magnetic field that is instantaneously producing a torque 
on that cell. The same procedure is applied to all cells, 
then, the demagnetization fields are recalculated, and the 
process is repeated iteratively until a desired tolerance is 
reached. A microscopic uniaxial anisotropy energy is also 
included, although using a strength that would be typi- 
cal for Permalloy, it is almost irrelevant when compared 
to the exchange and demagnetization effects. We have 
checked that this procedure gives the same final states as 
integration of the Landau-Gilbert equations with damp- 
ing. 

The internal magnetic energy E- lrA of the particle is 
calculated. This is the total magnetic energy minus 
the interaction energy with the applied magnetic field, 
—/I - H cxt - An applied magnetic field is used in the calcu- 
lations to move the net magnetic moment around, while 
it as well maps out the hysteresis loop. In one set of 
simulations, the hysteresis loop was calculated with the 
applied field axis within the xy-pleme at some angle 4>h 
to the x-axis. There, the magnetization makes an an- 
gle <p m to the x-axis. Then the internal energy could be 
found as a function -Ei n t(</>m), from which the anisotropy 
constant K\ is determined, by fitting to (U} in the form, 

£int (^m ) = Eq + K X sin 2 (f> m . (2) 

In another set of simulations, the applied field was set 
in the xz-plane, at some angle 9h to the x-axis. This 
tilts the net magnetic moment towards the z-axis by an 
angle 9 m from the cc-axis. Thus it gives E- m t(9 m ), which 
depends on both constants K\ and K3, according to 

E int (9 m ) =E + {K x + K 3 ) sin 2 9 m . (3) 

This allows anisotropy constant K3 to be determined 
from the net stiffness, K13 = K\ + K3. It is impor- 
tant to note, that these potential functions Ei nt (4>m) and 
E- m t(9 m ) found this way do not depend on the particular 
angle chosen between the applied field and the x-axis. 

In the following sections the Hamiltonian and algo- 
rithm is further specified. Some details about the demag- 
netization field calculation are given, especially concern- 
ing the Greens function. Finally the results for elliptic 
particles are discussed. 

II. THE PARTICLE MODEL AND ITS 
ENERGETICS 

We consider thin elliptical particles with dimensions 
L x x L y x L Zl where L x and L y refer to the major and 
minor diameters for the elliptical particles. The approach 
for a thin system has been presented in Rcf. 1 -; some of 
the main points towards finding the spatial structure of 
magnetization M(r) and the particle's internal energy 
are summarized here. 

The system is partitioned into cells of size a x a x L z 
on a square lattice grid, where there is saturation mag- 
netization M s within each cell. Thus a selected cell i 
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has a magnetic moment ir^ = M s a 2 L z rhi, that points in 
the direction of the unit vector rrij and has magnitude 
Mccii = M s a 2 L z . There is only a single layer of cells used, 
under the assumption that the perpendicular demagne- 
tization effect leads to magnetization nearly independent 
of z, for the thin systems under consideration. 

The exchange interaction in continuum theory is taken 
in terms of the exchange stiffness A (about 13 pJ/m for 
Py) as a volume integral, 



He 



A / dV Vm • Vm. 



(4) 



where rh — M/M s is the local reduced magnetization. 
When expanded on the square lattice of cells, this is 
equivalent to a nearest neighbor exchange term for the 
cells, 

H cx = —J^^ rhi ■ rhj, J = 2AL Z . (5) 

A uniaxial anisotropy energy K (about 100 J/m 3 for Py) 
is included as another volume integral 

■H uni = -K I dV (rh- uf -> -Ka 2 L z ^ (m, • u) 2 , (6) 

where the anisotropy axis here is taken as u = x. The 
externally applied magnetic field involves an energy of 
— B ■ fl for any dipole, so 

Hb = - J dV /Lto-ffext • M -> -p M s a 2 L z ^ H cxt ■ rh l . 

i 

Finally, the most important part of the interactions is the 
demagnetization field or dipolar interaction. Once the 
cells are defined on the grid with lattice spacing a, their 
dipole interaction could be described by a Hamiltonian, 
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[3(fhj ■ fij)(rhj ■ fjj) - <rhj ■ ™j] 



(8) 

However, this does not take into account the fact that 
the system is very thin. The demagnetization field can 
be found very accurately for thin systems using a Greens 
function approach^ 3 -. To do that, instead we start from 
the continuum dipolar energy, 



Wdd = -^o / >'Y 11, < ■ M 



(9) 



where Hm = — V$m is the demagnetization field at some 
point, and $m is its corresponding scalar magnetic po- 
tential. That field is produced by all the dipoles, accord- 
ing to a Poisson equation for magneto-statics: 



V 2 $m = Pm, where p M = —V • M. 



(10) 



Further, the discontinuity at the surfaces of the parti- 
cle can be modeled as a magnetic surface charge den- 
sity, <7m — M ■ h, where n is the outward normal. In 



particular, that gives charge densities of opposite signs, 
<jm = ±Mz on the upper and lower faces at z = 0, L z , 
respectively, under the assumption of uniform magneti- 
zation not depending on z within the cells. The field of 
those surface charges is responsible for keeping the mag- 
netization close to the xy plane. There are also surface 
magnetic charges at the edges of the island but those 
can be included into a localized volume charge for the 
cells there. But whether the magnetic charges are sur- 
face charges or volume charges makes no physical differ- 
ence, however, and the solution of the Poisson equation 
is formally 
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This can be used to find the demagnetization field at the 
point r = (x,y,z), and then averaging that result over 
z from z = to z = L z . The resulting demagnetization 
field at a cell centered at (x,y) has a vertical component 
Hm, z and some in-plane component HM,xy These are 
given by convolutions with appropriate 2D Green's func- 
tions, involving only the in-plane position, denoted here 
as f = (x,y). For the vertical demagnetization compo- 
nent, one gets 



H M , z (r) = / d A f' G z (r-r') M z (r'), r = (x,y) (12) 
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f 2 = x 2 + y 2 (13) 



Vr 2 + L 2 \r\_ 
For the in-plane components, there is 

H M , X y{?) = Id 2 ?' G xy (f-f' ) p(f). (14) 



G X y(f) 



2ttL 2 




(15) 



When applied, the unit vector points from source 

point f ' towards observation point f. There is a singu- 
larity in G z (r) as f — > 0, which is handled by averaging 
G z over a region with the same area as the cells being 
used, see Refj^ for further details on this averaging of 
the Green's functions. 

Together with appropriate finite-difference approxima- 
tions for the magnetic charge density, these expressions 
are used to get the demagnetization field. The actual 
evaluation of these convolution integrals was performed 
as multiplication in reciprocal space using a 2D fast 
Fourier transform approach 15 . To simulate a free par- 
ticle without the periodicity effects of the FFT (i.e., to 
avoid the wrap-around problem), the grid of the FFT is 
padded to a net size of N x x N y , where N x and N y are 
the smallest powers of two satisfying N x > 2L x /a and 
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N y > 2L x /a. Because we consider elongated particles, 
the calculations can be very fast due to N y being rather 
small compared to N x , for high aspect ratio particles. 
The FFT approach ends up giving the demagnetization 
field at the cell center positions (as well as at other points 
outside the particle, due to the padding). 

It is convenient to measure magnetic fields H ext and 
Hm in units of the saturation magnetization M s , just as 
done for the magnetization, rh = M/M s , and define the 
dimensionless fields, 



H 



The basic (and largest) unit of energy is the exchange 
J between neighboring cells. Then the total effective 
Hamiltonian can be written in units of J as 




(17) 
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FIG. 1: (Color online) The in-plane potential of an elliptical 
particle with a 5:1 aspect ratio, mapped out while determining 
the hysteresis loop (Fig. [3}. The angle <j> m is the direction of 
the net particle moment jl in the easy-plane. The points come 
from simulations at the different indicated angles 4>h of the 
applied field from the long axis; all fall onto the same curve. 
The fit gives a reliable estimate of anisotropy constant K\. 



This is written 
change length A, 
anisotropy k, defined as 



in terms of the ferromagnetic ex- 
,. x and the scaled dimensionless uniaxial 
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The magnetic internal energy E{ nt is of most interest. 
That is this Hamiltonian, but with the interaction with 
the external magnetic field (the last term) removed. 

For the calculations we used the values for Permalloy, 
M s « 860 kA/m, A ps 13 pJ/m, K rj 100 J/m 3 , then 
these give A ex ~ 5.3 nm and K « 1.1 X 10~ 4 . Due to this 
small value of n, in most of the calculations the intrinsic 
uniaxial anisotropy energy is negligible compared to the 
other energies of the system. In most of the simulations 
the cell size was a = 2.0 nm, except for the smallest high 
aspect ratio particles, where values as low as a = 0.5 nm 
were used, to produce a smoother edge to the particle. 
These are sufficiently less than the exchange length to 
give a reliable description of the internal magnetic struc- 
ture. 



That local field enters the undamped dynamic equation 
of motion, 



drfii 



jrhi x Bi, (19) 
and it is given by the variation of the Hamiltonian, 
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Alternatively, this is the same as 
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(21) 



nbrs 



+ 



2n{rhi ■ u)u + -h M ,t + h ext 



Either way, these define a unit vector along which to 
point the magnetic moment of cell i, rhi — > rh';, where 



III. CALCULATION PROCEDURES 

The iteration procedure that moves the system towards 
the nearest local equilibrium is a local spin relaxation 
algorithm 1 -, that points each cell's magnetic moment rhi 
(or its unit vector rhi) along its local magnetic field Bi. 



rh 1 ; = bi 



_B_ 

\Bi 



(22) 



The alignment of rhi parallel to bi is performed for every 
site of the grid, after which the demagnetization field 
must be recalculated. The process moves the system 
towards lower energy. Each cell would stop moving if 
all became simultaneously aligned self-consistently with 
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FIG. 2: (Color online) The out-of-plane potential of the same 
elliptical particle with a 5:1 aspect ratio as in Fig. [T] mapped 
out while determining the hysteresis loop (Fig. [4]). The angle 
9 m is the tilting of the net particle moment fl out of the 
easy-plane. The points from simulations at different angles 
6h of the applied field are combined into one curve. The fit 
gives a reliable estimate of the combined anisotropy constant 
K13 — Ki + K3. 

their local magnetic fields. This does not insure a global 
energy minimum, however, and the procedure does have 
the capability to generate the hysteresis loops. The iter- 
ation is started from a partially aligned state of the cell 
dipoles, which are given some small random fluctuations 
away from perfect alignment. For the hysteresis calcula- 
tion, though, the last relaxed state at one applied field is 
the initial state for the next value of applied field. 



As the iteration proceeds, periodically (every 200 iter- 
ations of the system) the total magnetic moment fl of the 
particle is calculated, by summing over the cell dipoles, 

pL=(n x ,Hy, n z ) = /i C eii y] rhj (23) 

i 

The iteration is stopped when the changes in any com- 
ponent of fx are less than 1 part in 5 x 10 7 for two states 
separated by 200 iterations. This is actually a more strin- 
gent stopping requirement than waiting for the energy to 
converge to the same precision. 

In one set of simulations, the applied magnetic field 
was directed within the xy-plane at an angle 4>h to the x- 
axis [-ff cx t = -ff e xt(cos 4>h, sin <pjj , 0)]. This results in the 
net magnetic moment /2 lying within the xy-plane, due to 
the strong perpendicular demagnetization. The moment 
makes some angle <f> m < <Ph to the x-axis, calculated 
from <p m = tan _1 (/Ltj / //x x ). Thus, the in-plane potential 
energy function E- lnt (4> m ) could be found, by scanning 
in applied magnetic field strength along the chosen axis, 



and thereby calculating the hysteresis loop. Then K\ was 
found by fitting to the form in 

In the other set of simulations, the magnetic field was 
applied tilting out of the xy-plane, making an angle 9h to 
the cc-axis, that is, H ext — iJ cxt (cosf?y, 0, sin##). This 
pulls ft up an angle 9 rn < 8 H from the easy (xy) plane, 
where 6 m = taxi (jj, z / fi x ) , and gives the opportunity to 
measure the potential E- lni {9 m ). Again, the internal en- 
ergy is calculated from the total Hamiltonian minus the 
interaction term with the applied external field. That en- 
ergy was fitted to the out-of-plane potential ([3]), whose 
stiffness is due to the combination, K\% — K\ +K3. Fur- 
ther, the potentials obtained did not depend on the choice 
of 4>h or Oh- This could be seen by combining the inter- 
nal energy curves for applied field at 45° and 80° to the 
cc-axis. 

For the hysteresis curves, the total magnetic moment 
\x was calculated, and normalized by the particle volume 
V to get the averaged magnetization inside the particle, 
(M) = Then the component of (M) along the applied 
field axis is found, 

(M h ) = (M) ■ h cxt (24) 

After scaling by the saturation magnetization, this is 
plotted versus the applied field magnitude also scaled by 
saturation magnetization (h cxt — H cxt /M s ). 



IV. RESULTS FOR ELLIPTICAL PARTICLES 

We considered thin elliptical particles with thicknesses 
all l/20 th of the length, i.e., 33 = L x /L z = 20, and 
aspect ratios g\ = L x / L y — 3, 5, and 8. The lengths 
ranged from 120 nm to 480 nm. Some typical results for 
the internal energy curves are shown in Fig. [T]for the in- 
plane potential and Fig. [5] for the out-of-plane potential 
of an elliptical particle with g\ — 5, with major axis 
240 nm, minor axis 48 nm and thickness 12 nm. The 
potentials for in-plane motion of /2 fit very well to the 
functional form, 

E int ((f> m ) =E a + K 1 sin 2 cj> m . (25) 

The constant Eq is an irrelevant ground state energy 
when the particle is magnetized along its long axis. This 
same form also applies to the potential E- m t(Q m ), but 
with a coefficient K\% = K\ + K3. The fits are best for 
smaller particles, where the cells stay strongly aligned 
with each other, and the reversal can be considered close 
to a uniform rotation process, for the most part. For 
the larger particles (length > 400 nm) this global align- 
ment is lesser and the fits are good but with considerably 
greater \ 2 - Even so, the internal magnetization structure 
of the relaxed states tends to be close to uniform. 
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Li = 120 nm 


L x = 240 nm 


L x = 480 nm 


gi = 2 


ifi = 6.35J 


27.3J 


111J 


gi = 3 


7.32J 


31.9J 


134J 


ffi = 5 


6.96J 


31.5J 


133J 


3i = 8 


7.39J 


29.5J 


118J 



TABLE I: Values of the in-plane anisotropy constant K\ in 
units of J = 2AL Z for the indicated particle sizes and aspect 
ratios gi — L x /L y . All of the particles calculated have gz = 
L x /L z = 20. 





La: = 120 nm 


L x = 240 nm 


L x = 480 nm 


gi = 2 


#13 = 79.0J 


314J 


1250J 


gi = 3 


50.7J 


201J 


804J 


ffi = 5 


28.1J 


111J 


444J 


9i = 8 


15.7J 


62.6J 


250J 



TABLE II: Values of the combined anisotropy constant Kiz 
in units of J = 2AL Z for the indicated particle sizes and 
aspect ratios gi = L x /L y . All of the particles calculated have 
<?3 = L x /L z = 20. 



Fitting results are summarized in Tables HI and [TT1 with 
the constants presented in units of J = 2 AL Z . The hard- 
axis anisotropy constant K3 derived from those results is 
presented in Table Mil The constant K3 is consistently 
stronger than the easy-axis constant Ki, as to be ex- 
pected from the greater surface area of the lower and 
upper faces at z = 0,L Z , compared to the very limited 
surface area of the edge of the ellipse. The energy unit 
J itself varies according to the thickness. Thus it makes 
sense to also look at results for the constants in joules. 
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FIG. 3: (Color online) Hysteresis loops for an elliptical parti- 
cle as in Fig. [T]with an in-plane applied field at the indicated 
angles cj>n to the long axis of the particle. 




FIG. 4: (Color online) Hysteresis loops for an elliptical parti- 
cle as in Fig. f2]with the applied field tilted out of the sy-plane 
at the indicated angles 9h from the long axis of the particle. 





L x = 120 nm 


L x = 240 nm 


L x = 480 nm 


gi = 2 


K A = 72.7 J 


287J 


1140J 


gi = 3 


43.4J 


169J 


670J 


gi =5 


21. 1J 


79.9J 


311J 


gi = 8 


8.30J 


33.1J 


132J 



TABLE III: Values of the hard-axis anisotropy constant K3 
in units of J = 2AL Z for the indicated particle sizes and 
aspect ratios g\ — L x /L y . All of the particles calculated have 
gz = L x /L z = 20. 



Generally, K3/J increases proportional to the area of 
the ellipse, jTtL x L v , multiplied by the thickness l z) so 
that in fact K3 (in joules) is linearly proportional to 
the volume of the particles. Also, one sees that K3 de- 
creases with increasing aspect ratio for particles of the 
same length; this is because the particle volume is de- 
creasing. On the other hand, K\ / J depends very weakly 
on the aspect ratio for the particle sizes tested. In addi- 
tion, the calculations suggest that K\ increases somewhat 
faster than the particle volume. The weak dependence of 
A'i on the shape of the ellipse (at these larger values of 
gi) is surprising. 

To clarify the results we also show the constants 
converted to energy densities, both K\JV and K3/V 
in joules/nm 3 , in Fig. [SJ The actual units are the 
exchange stiffness A (units of joules/nm) divided by 
squared nanometers. One finds very little dependence 
of either energy density constant, K/V, on the particle 
size, however, again it is clear that K3 is always larger 
than K\. Furthermore, the easy-axis anisotropy constant 
Ki/V does increase rapidly with the in-plane aspect ratio 
gi, and the relation could be close to a linear relation- 
ship. Although the values of K3/V are always greater 



7 



0.03 
0.025 

— 0.02 


a 

< 0.015 
> 

* 0.01 
0.005 



1 1 1 1 1 1 
K./V (dashed) 


i i I 
g,=2 


I i i | i i i i 

All for g=20 












- 






















-<J? 




▲ 




A. J 


K/V (solid) 


^ 

,f,=5 




- ♦ 


_g,= 3 













Too 



200 



300 
L (nra) 



400 



500 



FIG. 5: (Color online) The anisotropy constants K\ (solid 
curves) and K3 (dashed curves) scaled by elliptical particle 
volume, versus particle lengths, for the indicated gi aspect 
ratios. All data has g% = 20. The values of K/V are given 
in units of A/nm 2 , where A is the exchange stiffness. K\/V 
increases with aspect ratio while K3/V decreases, and they 
become equal at high aspect ratio. 



than the corresponding Ki/V, these hard-axis energy 
densities K 3 /V decrease slightly with increasing aspect 
ratio g%. At large aspect ratio, the two constants become 
nearly the same, which would have to be the case for a 
needle-shaped magnet. 



A. The magnetization structure 

In the high aspect ratio particles, the magnetization 
states are very close to uniform, even when undergoing 
the reversal. The elongated particle has such a strong 
anisotropic effect that the magnetization cells move al- 
most in a synchronized motion. For particles with smaller 
aspect ratio, one starts to see some weak variations in the 
magnetization inside the particle. 

To get an idea of the size of this effect, some configu- 
rations are presented for ellipses with g\ = 2, which has 
the strongest effect of all the particle shapes presented 
earlier. In Fig. |6] some configurations are shown for a 
120 x 60 x 6 nm particle, at different applied field strengths 
45° to the particle's long (+x) axis. The points shown 
are at (a) close to saturation, (b) zero applied field, (c) 
a negative field close to reversal, and (d) a negative field 
just after reversal. For the most part, the magnetization 
stays nearly uniform for this relatively small particle. 

Another example is presented in Fig. [7J like the first 
example, but 2x larger in all three dimensions. The 
four configurations shown correspond to the same four 
types of states as presented for the smaller particle. The 
main difference here is that a nonuniform magnetic struc- 
ture develops. At zero field, the structure points in- 
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FIG. 6: Magnetic configurations for a 120 x 60 x 6 nm par- 
ticle with magnetic field applied at +45° above a horizontal 
axis pointing to the right. The arrows are the coarse-grained 
averages of 3 x 3 groups of cells. In (a), the external field is 
h = 0.20; in (b) h = 0.0; (c) h = -0.030, just before reversal; 
(d) h — —0.032, just after reversal. 



ward/outward towards the poles on the long axis. For 
the configurations just before and after reversal, a wave- 
like structure is present. These spatial variations are due 
to the dipolar interactions; in even lower aspect ratio par- 
ticles (gi < 2), they lead to C-states and even vortices 
entering the particle. 
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FIG. 7: Magnetic configurations for a 480 x 240 x 24 nm par- 
ticle with magnetic field applied at +45° above a horizontal 
axis pointing to the right. The arrows are the coarse-grained 
averages of 9 x 9 groups of cells. In (a), the external field is 
h = 0.20; in (b) h = 0.0; (c) h = -0.025, just before reversal; 
(d) h — —0.027, just after reversal. Note the enhanced cur- 
vature of the field compared to that in the smaller particle in 

Fig. m 



B. Particles with lower aspect ratio gi < 2 

When gi — > 1, the ellipse becomes circular and the 
easy-axis anisotropy must vanish. Using smaller g\ 
is a way to produce particles with weaker easy-axis 
anisotropy constant. However, as the system becomes 
closer to circular, the lowest energy configuration, es- 
pecially near zero applied magnetic field, tends to be 
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FIG. 8: (Color online) The anisotropy constants K\ (solid 
curves) and K3 (dashed curves) scaled by elliptical particle 
volume, versus particle thicknesses, for the indicated g\ aspect 
ratios. All the data is for particles of length L x = 240 nm. 
The K3/V constant crosses below zero for the thickest high- 
aspect-ratio particles, which have become needle-like and no 
longer satisfy the assumption of a thin particle. That is the 
case of a particle with only uniaxial anisotropy. 



nonuniform. The ground state can tend towards a C- 
state or a vortex state if the particle is of sufficient size. 
The above results do not apply to that situation, espe- 
cially because the nonuniform magnetization cannot be 
mapped into the model of an individual magnetic mo- 
ment moving in an effective potential. 

To verify this, some particles were also calculated at 
small cllipticity, where K\ 0, using g\ = 1.25 and 
<7l = 1.11. Generally, at these ratios, if there was a stable 
single-domain ground state (for smaller particles only), 
the tendency is for the moments to try to follow the bor- 
der, and point inwards or outwards from the poles at the 
long ends. At larger particle size this tilting eventually 
moves the system irreversibly to a vortex ground state. 
Until the vortex state is reached, an effective potential 
can be estimated, however, from the practical point of 
view it may be of limited use. 



C. Thicker particles 

The particles with 53 = 20 can be too thin to hold a 
magnetic moment stable against room-temperature ther- 
mal fluctuations. Thus it is important to consider the 
changes when thicker particles are used. Further cal- 
culations were carried out for 240 nm long particles to 
get results for g 3 = 20, 15, 10, and 8, corresponding to 
thicknesses of 12, 16, 24 and 30 nm, respectively. The 
results for K\/V and K3/V are shown in Fig. HI As 
could be expected, the thicker particles have weaker out- 
of-plane anisotropy K3/V, while Ki/V increases due to 
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the thicker lateral edges, but at a rate less than linear in 
the thickness. We expect that these per- volume energy 
constants have only very weak dependence on the parti- 
cle length, as was already seen in the results presented 
above for 12.0 nm thickness. 



V. CONCLUSIONS AND DISCUSSION 

The anisotropy properties of thin elliptical ferromag- 
netic particles have been estimated, based on a 2D mi- 
cromagnetics model that employs Green's functions for 
the calculation of the demagnetization fields. For the 
high-aspect-ratio particles being considered, the magne- 
tization was found to be close to uniform inside the parti- 
cles. Then it was possible to map out the changes in the 
internal energy versus the direction of the net magnetic 
moment ft, which itself acts as a collective coordinate. 
The typical particles tend to have stronger anisotropy 
in the hard-axis direction (K3/V) than in the easy-axis 
direction (K\/V), however, these two energy scales ap- 
proach each other for needle-like particles, as expected. 
The results could be of practical application in the design 
and analysis of artificial spin-ice with desired dynamics, 
beyond the usual Ising energetics. 

In the theoretical study of artificial spin ice materials, 
it is usual to replace the islands by point-like dipoles with 
an Ising-like behavior. Indeed, all theoretical calculations 
for the properties of these systems were obtained with 
this approach. However, a more realistic description of 
these artificial spin ices should require models beyond the 



Ising approximation, such as continuous magnetic mo- 
ments with anisotropy considered in this work. In such 
a case, although the main properties of a spin ice system 
may not undergo strong alterations, several quantities 
would change their values. For instance, a recent work 
about the thermodynamics of the square lattice^ has 
suggested a possible phase transition in this system, oc- 
curring at a temperature of 7. 2D, where D is the coupling 
constant of the dipolar interaction among the islands. Of 
course, the transition temperature or similar quantities 
should be dependent on the island sizes and anisotropics, 
but this dependence cannot be perceived with the Ising 
approach. It is very probable that the correct critical 
temperature must be much smaller than 7. 2D since the 
total magnetic moment of an island has more degrees 
of freedom, and effectively moves in a softer potential. 
In addition, the properties must also be dependent on 
the islands' shapes, etc. So, the results obtained here 
are of fundamental importance for developing this field 
not only theoretically but also experimentally, suggest- 
ing protocols for improving experiments, and including 
studies about their dynamics. 
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